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Abstract 

The  structure,  ionic  conductivity  and  strength  of  poly(ethylene  oxide)  (PEO)  sulfonic  acid-based  solid  polymer  electrolytes  with  various 
contents  of  the  PEO  polymer  are  analyzed  using  molecular  dynamics  simulations.  To  quantify  the  electrolyte  structure,  comprehensive 
coordination  and  dimensional  analyses  are  carried  out.  Ionic  conductivity  is  determined  by  computing  the  Einstein-diffusion  based  conductivity 
for  all  ionic  particles  in  the  system.  The  strength  of  the  electrolyte  is  quantified  by  carrying  out  a  set  of  molecular  simulations  of  the 
uniaxial  deformation  process  under  constant  stress-rate  conditions.  The  results  obtained  indicate  that  increased  hydration  improves  the  ionic 
conductivity,  but  degrades  the  strength  in  solid  polymer  electrolytes  at  hand.  In  addition,  it  appears  that  there  is  an  optimal  level  of  the  PEO 
polymer  content  in  these  materials,  which  at  the  same  hydration  level,  yields  the  best  combination  of  ionic  conductivity  and  strength. 
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1.  Introduction 

Over  the  last  two  decades,  solid  polymer  electrolytes  have 
been  the  subject  of  intense  investigation  due  to  their  wide- 
range  use  as  non-corrosive  solid  electrolytes  in  various  elec¬ 
trochemical  applications,  electronic  equipments,  medical  de¬ 
vices  and  in  fuel  cells  for  electrical  and  hybrid  vehicles  [1-5] . 
A  comprehensive  review  of  the  synthesis,  properties  and 
applications  of  the  proton-conducting  solid  polymer  elec¬ 
trolytes  can  be  found  in  Refs.  [6,7]. 

In  a  comprehensive  study,  Herranen  et  al.  [8]  car¬ 
ried  out  synthesis  and  testing  of  a  promising  class  of 
solid  proton-conducting  polymer  electrolytes  based  on  the 
poly  (ethylene  oxide)  (PEO)  with  sulfonic  acid  end  groups. 
These  oligomeric  electrolytes,  referred  to  as  the  PEO  sul¬ 
fonic  acids,  are  embedded  in  a  high  molecular  weight  PEO 
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matrix  to  obtain  mechanically  stable  electrolyte  membranes 
with  ionic  conductivities  as  high  as  1.5  mS/cm.  The  molecu¬ 
lar  structure  and  the  interactions  between  the  molecules  and 
the  charge  carriers  were  characterized  using  thermal  anal¬ 
ysis,  impedance  spectroscopy,  NMR  and  IR  spectroscopy 
[8]. 

The  development  of  new  materials  via  purely  experimen¬ 
tal  means  is  a  time-consuming  and  costly  proposition.  This 
is  particularly  true  in  the  case  of  solid  polymer  electrolytes, 
since  the  mechanisms  controlling  various  properties  of  these 
materials  (primarily  ionic  conductivity,  hydration  behavior, 
osmosis,  etc.)  are  not  well  understood  and,  hence,  offer  very 
little  guidance  to  the  material  design  process.  Fortunately,  a 
rapid  increase  in  the  computational  resources,  and  the  devel¬ 
opment  of  new  highly  efficient  software  packages  combined 
with  parallel  advances  in  the  experimental  materials  char¬ 
acterization  tools  now  offer  new  opportunities  for  obtaining 
an  insight  into  the  atomic-level  structure  and  mechanisms  in 
these  materials. 
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At  present,  there  is  a  great  need  for  solid  polymeric  elec¬ 
trolytes  with  an  ionic  conductivity  in  a  mS/cm  range,  at  tem¬ 
peratures  around  100  °C.  Among  such  solid  electrolytes  those 
which  are  proton-conducting  are  of  particular  importance  due 
to  their  potential  use  as  membranes  in  the  polymer  electrolyte 
membrane  (PEM)  fuel  cells.  Despite  a  significant  research, 
which  was  devoted  to  understanding  conductivity  in  proton¬ 
conducting  polymers,  the  underlying  mechanism  of  this  pro¬ 
cess  is  still  not  well  understood.  It  is  generally  accepted  that 
there  are  two  major  contributions  to  the  proton-conductivity 
in  solid  polymeric  electrolytes  in  the  presence  of  water:  (a) 
diffusion  of  individual  protons  (or  more  precisely  hydronium 
ions,  H30+)  through  the  polymer  electrolyte  medium  and 
(b)  a  proton  hopping  process  within  which  protons  combine 
with  the  water  molecules  to  form  hydronium  ions  and  sub¬ 
sequently  hop  from  these  hydronium  ions  to  the  surrounding 
water  molecules  [9] .  While  the  proton/hydronium  diffusion 
mechanism  can  be  studied  using  classical  molecular  dynam¬ 
ics  simulations,  and  the  corresponding  value  of  the  diffusion 
coefficient  calculated  from  the  simulation  results,  an  explicit 
treatment  of  the  proton  hopping  process  entails  the  use  of 
a  quantum  mechanical  approach.  Consequently,  the  proton 
hopping  mechanism  can  be  studied  only  in  relatively  small 
computational  systems  consisting  of  few  tens  of  particles  and, 
hence,  the  contribution  of  proton  hopping  to  the  overall  pro¬ 
ton  transfer  process  in  macroscopic  systems  may  not  be  prop¬ 
erly  assessed. 

In  a  series  of  papers,  Ennari  et  al.  [10-17]  proposed  a 
new  approach  to  dealing  with  the  proton  hopping  mecha¬ 
nism.  Within  this  approach,  a  new  particle  named  proton  has 
been  created  and  assigned  the  mass  and  the  charge  of  a  stan¬ 
dard  proton.  However,  the  van  der  Waals  radius  of  this  par¬ 
ticle  has  been  derived  from  the  condensed-phase  properties 
of  polymers  containing  highly  polarized  hydrogen  (e.g.  wa¬ 
ter  or  methanol).  Consequently,  the  new  particle  can  react 
with  the  surrounding  molecules  via  strong  Coulomb  elec¬ 
trostatic  forces  and,  hence,  jump  from  one  water  molecule 
to  the  other.  In  this  way,  molecular  dynamics  simulations 
can  be  used  to  assess  the  contribution  of  the  proton  hopping 
mechanism  to  the  proton  transport  in  solid  polymeric  elec¬ 
trolytes.  However,  it  should  be  noted  that  while  the  approach 
proposed  by  Ennari  et  al.  [10-17]  is  capable  of  mimicking 
the  proton  hopping  mechanism,  it  does  not  represent  the  true 
proton  hopping  process  in  which  chemical  bonds  are  broken 
within  the  hydronium  molecules  and  generated  within  the 
water  molecules.  Nevertheless,  the  approach  of  Ennari  et  al. 
[10-17]  offers  a  good  first  step  toward  dealing  with  a  very 
complex  problem,  and  since  it  yields  diffusivity  values  com¬ 
parable  to  their  experimental  parts  [10-17],  it  is  adopted  in 
the  present  work. 

The  objective  of  the  present  work  is  to  extend  the  ap¬ 
proach  of  Ennari  et  al.  [10-17]  in  order  to  explore  the  trade¬ 
off  between  a  hydration-induced  increase  in  ionic  conductiv¬ 
ity  in  the  PEO  sulfonic  acid-based  solid  polymer  electrolytes 
and  the  accompanied  decrease  in  the  material  strength.  When 
these  solid  polymer  electrolytes  are  used  in  electrochemical 


membrane  applications,  they  must  generally  meet  require¬ 
ments  pertaining  to  both  the  levels  of  their  ionic  conductivity 
and  their  strength. 

The  organization  of  the  paper  is  as  follows:  details  of  the 
computer  simulation  procedure  and  the  methods  employed  to 
analyze  materials  microstructure  and  properties  of  the  mate¬ 
rial  studied  in  the  present  work  are  overviewed  in  Section  2. 
The  main  results  obtained  in  the  present  work  are  presented 
and  discussed  in  Section  3 .  The  key  conclusions  resulted  from 
the  present  study  are  summarized  in  Section  4. 

2.  Computational  procedure 

2.7.  Constituents 

In  all  the  calculations  carried  out  in  the  present  work,  the 
following  constituents  are  used:  (a)  PEO  molecules  based 
on  a  (CH2CH2O)  repeat  unit;  (b)  PEO  sulfonic  acid  di¬ 
anions  with  a  formula  “03SCH2CH2CH2  0(CH2CH20)4 
CH2CH2CH2SO3-  and  a  molecular  weight,  Mw  =  436;  (c) 
water  molecules  (H2O);  (d)  hydronium  ions  (H30+);  (e)  pro¬ 
tons  (H+).  A  schematic  of  these  constituents  and  the  charges 
of  their  atoms  are  displayed  in  Fig.  1. 

2.2.  Construction  of  the  computational  cells 

To  explore  the  effect  of  the  PEO  polymer  content  on  ionic- 
transport  and  mechanical  properties  of  PEO  sulfonic  acid- 
based  solid  electrolytes  at  the  same  weight-percent  level  of 
water,  three  (cubic)  cells  are  constructed. 

Each  of  the  three  cells  contains  four  PEO  sulfonic  acid 
dianions  with  a  formula,  “03SCH2CH2CH2  0(CH2CH20)4 
CH2CH2CH2S03“  (molecular  weight,  Mw  =  436),  four  (H+) 
protons  and  four  (H30+)  hydronium  ions.  To  explore  the  ef¬ 
fect  of  the  PEO  polymer  content,  the  three  cells  contained  one 
PEO  molecule  with  the  degree  of  polymerization:  of  40  (Cell 
A),  60  (Cell  B)  and  80  (Cell  C).  To  attain  a  fixed  40wt.% 
of  water,  the  three  cells  are  also  assigned  80,  111  and  142 
water  molecules,  respectively.  The  edge  lengths  of  the  three 
unit  cells  are  next  adjusted  to  minimize  the  potential  energy 
in  each  case.  This  procedure  resulted  in  the  following  mate¬ 
rial  densities:  0.90,  0.96  and  1.04  g/cm3  for  the  three  cells, 
respectively.  Furthermore,  for  each  type  of  cells,  50  configu¬ 
rations  are  constructed. 

In  the  present  work,  the  structure  and  properties  of  PEO 
sulfonic  acid-based  electrolytes  is  studied  using  a  larger 
number  of  initial  configurations  and  relatively  short  runs 
(1000  ps).  This  is  done,  because  it  was  observed  that  the  start¬ 
ing  structures  generated  by  the  Amorphous  Cell  Builder  do 
not  change  substantially  during  longer  runs.  Hence,  better 
statistics  is  expected  from  a  larger  number  of  short  runs  than 
from  fewer  long  runs.  In  each  case,  the  final  properties  are  cal¬ 
culated  using  only  the  trajectory  results  from  the  last  200  ps. 

Individual  constituent  molecules  (e.g.  PEO,  H30+,  etc.) 
are  constructed  using  the  Accelrys’  software  package 
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Fig.  1.  Partial  changes  for  the  atoms  in:  (a)  PEO;  (b)  PEO  sulfonic  acid  dianion;  (c)  water;  (d)  hydronium;  (e)  proton. 


Materials  Visualizer  [18],  while  the  three  cells  containing 
the  constituents  described  above  are  constructed  using  the 
Accelrys’  package  Amorphous  Cell  Builder  [19].  In  all  sim¬ 
ulations,  the  periodic  boundary  conditions  are  applied  in  all 
three  directions  perpendicular  to  the  faces  of  the  cells. 

2.3.  Forcefields 

All  the  calculations  carried  out  in  the  present  work  are 
done  using  the  PCFF  forcefield  [20]  and  its  recent  modifi¬ 
cation  called  NJPCFF  forcefield  [21]  that  enable  computa¬ 
tion  of  the  interactions  between  ions  and  atoms/molecules. 
The  PCFF  forcefield  contains  11  intra-molecular  bonding 
terms  such  as  a  single-bond  stretching  term,  a  two-bond  angle 
term,  a  dihedral  middle-bond  torsion  term,  etc.  In  addition, 
the  forcefield  includes  two  intra-  and  inter-molecular  non¬ 
bonding  (van  der  Waals  and  Coulomb)  terms.  Within  the 
NJPCFF  forcefield,  S— O-  bond  stretching  and  “O— S-O- 
two-bond  angle  parameters  are  defined  for  the  PEO  sulfonic 
acid  dianion.  In  addition,  a  proton  particle  is  introduced  with 
a  mass  of  1.00797w  ( u  is  the  atomic  mass  unit),  a  charge  q 


normalized  by  the  electronic  charge  e  of  1.0,  and  the  van  der 
Waals  non-bonding  parameters  r*  =  2.5  and  Si  =  0.013. 

2.4.  Selection  of  the  initial  configurations 

For  each  of  the  150  configurations  of  the  three  types 
of  cells  described  in  Section  2.2,  the  energy  is  initially 
minimized  for  10,000  steps  using  the  steepest-decent 
molecular-mechanics  method.  During  these  and  the  subse¬ 
quent  calculations,  the  Cell  Multipole  summation  method 
with  a  relative  dielectric  constant  of  1.0  and  an  update  width 
of  1  A  is  used  when  dealing  with  the  non-bonding  (van 
der  Waals  and  Coulomb)  interactions.  Ten  configurations 
for  each  of  the  three  cells  with  the  lowest  energies  are 
selected  for  the  subsequent  analysis.  The  energy  of  the 
selected  cells  is  next  minimized  using  a  combination  of 
the  steepest-decent  and  the  conjugate-gradient  methods 
until  a  convergence  defined  by  the  maximum  molar  energy 
derivative  of  10  kcal/(mole  A)  is  reached. 

The  resulting  atomic  configurations  are  next  subjected 
to  the  standard  molecular  dynamics  simulations  for  1000  ps 
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Fig.  2.  An  example  of  the  computation  cell  containing  the  PEO  polymer 
with  the  polymerization  index  of  40  in  the:  (a)  original  configuration;  and 
(b)  in  the  deformed  state. 

with  a  fixed  time  increment  of  1  fs.  The  configurations  are 
treated  as  NVT  ensembles  within  which  the  temperature  is 
held  at  298  K  using  the  Anderson  thermostat  method  [22]. 
The  Cell  Multipole  summation  method  with  a  relative  di¬ 
electric  constant  of  1.0  and  an  update  width  of  1  A  is  again 
used  to  efficiently  handle  the  non-bonding  interactions. 

The  molecular  statics  (energy  minimization)  and  the 
molecular  dynamics  simulations  described  above  are  carried 
using  the  Accelrys’  software  package  discover  [23].  An  ex¬ 
ample  of  the  initial  cell  of  type  A  is  given  in  Fig.  2(a). 

2.5.  Vibrational  spectral  analysis 

A  comparison  of  the  computed  and  measured  vibrational 
spectra  is  generally  used  to  validate  the  forcefield  used  in 
molecular  modeling.  The  vibrational  spectrum  of  a  molecule 
can  be  calculated  using  the  local  mode  method.  Within  this 
method,  low-velocity  components  of  the  constituent  parti¬ 
cles  are  calculated  using  the  classical  molecular  dynamics 
simulations.  To  calculate  the  fast  vibrations  of  the  atoms,  on 
the  other  hand,  several  representative  “equilibrium”  atomic 
configurations  are  first  generated  using  molecular  dynamics 
simulations.  For  each  of  these  configurations,  the  molecular 
vibrational  frequencies  are  then  determined  using  quantum 
mechanical  calculations  within  which  each  particle  is  treated 


as  an  quantum  oscillator  embedded  in  an  effective  medium 
of  the  remaining  atoms.  The  local  potential  energy  function 
for  such  an  effective  medium,  Veff ,  is  defined  as  the  following 
polynomial: 

Vefi(Q)  =  K0  +  K2Q2  +  K3Q3  +  K4Q4  (1) 


where  Q  defines  coordinates  of  the  species  resulting  from  the 
molecular  dynamics  simulation  runs. 

Next,  the  following  one-dimensional  Schrodinger  equa¬ 
tion  is  solved: 


n 

2 M 


d2 

d Q2 


+  veff(G) 


V(Q)  =  E9(Q) 


(2) 


where  M  is  the  standard  (molecular  geometry)  G  matrix  in  a 
normal  mode  analysis  [26]  and  ^(Q)  is  the  wave  function. 
The  accuracy  of  the  frequencies  calculated  using  the  local 
mode  method  is  generally  considered  to  be  around  100  cm-1 . 


2.6.  Conformational  and  structural  analyses 

To  characterize  the  intra-molecular  and  the  inter- 
molecular  structures  of  the  materials  studied  in  the  present 
work,  the  following  quantities  are  computed  and  analyzed:  (a) 
pair  correlation  functions;  (b)  coordination  numbers;  (c)  radii 
of  gyration;  (d)  end-to-end  molecular  distances;  (e)  charac¬ 
teristic  molecular  ratio. 


2.6.1.  Pair  correlation  functions 

Pair  correlation  functions  (often  referred  to  as  radial  dis¬ 
tribution  functions),  gXm .  mZ(r),  define  the  probability  of  finding 
a  pair  of  atoms  (labeled  as  x  and  z)  at  distance  r  apart,  rela¬ 
tive  to  the  corresponding  probability  in  the  same  material  (at 
the  same  density),  but  under  the  condition  of  a  completely 
random  distribution  of  the  particles  [24] . 


2.6.2.  Coordination  number 

The  coordination  number  for  particle  z,  nXm .  mZ(r),  can  be 
calculated  from  the  corresponding  pair  correlation  functions 
as: 

nx...z(r)  =  47T777T  [  8x..As)s2ds  (3) 

where  nXm .  .z(r)  represents  the  number  of  x  particles  located 
within  a  sphere  centered  at  the  particle  z  and  having  a  radius 
r,  Nz  is  the  total  number  of  z  particles  in  the  system,  while 
(V)  is  the  average  system  volume. 

2.6.3.  Radius  of  gyration 

The  size  of  a  polymer  chain  can  be  represented  using  its 
radius  of  gyration  whose  radius,  R2g,  is  defined  as  an  average 
square  distance  of  the  constituent  particles  of  the  molecule 
from  the  molecule’s  center  of  gravity.  The  radius  of  gyration 
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can  hence  be  expressed  as: 

1  N 

R2g=NH«R‘-RG?)  (4) 

i=  1 

where  N is  the  number  of  constituent  particles,  Rj  the  position 
vector  of  particle  i,  ( )  denotes  an  average  and  the  center  of 
gravity  position  vector  is  defined  as: 

*G  =  ^X>  (5) 

i=  1 


the  transported  species  from  the  opposite  surface  of  the  mem¬ 
brane.  However,  it  is  well  established  that  transport  through 
the  membrane  is  the  slowest  and,  hence,  the  rate  controlling 
process  [25]. 

Within  an  atomistic-modeling  framework,  the  diffusion 
coefficient  for  species  a,  Da ,  can  be  defined  by  the  following 
equation  [24]: 


1  d  Na 

A*  =  ztt  lim  t  7 -  Ri 

6Na  t  y  oo  d t  ^ 
i—l 


(0)]2> 


(9) 


A  simple  manipulation  of  Eqs.  (4)  and  (5)  yields: 

i=  1  j=n+ 1 

2.6.4.  The  mean  square  end-to-end  molecular  distance 
An  alternative  measure  of  the  size  of  a  polymer  chain 
is  the  mean  square  end-to-end  distance,  ( R 2).  This  mea¬ 
sure  of  the  molecule  size  (as  well  as  the  characteristic  ratio, 
Cn)  takes  into  account  the  short-range  interactions  (as  mani¬ 
fested,  for  instance,  by  the  rotational  isomericity,  the  torsional 
chain  stiffness,  etc.)  between  various  segments  of  a  molec¬ 
ular  chain.  The  mean  square  end-to-end  distance  is  defined 
as: 

(R2)  =  (Rm-Ri)2  (7) 

where  the  subscripts  M  and  1  are  used  to  denote  the  last  and 
the  first  particle  in  the  backbone  of  a  molecular  chain. 


2.6.5.  Characteristic  ratio 

The  chain  dimensions  can  also  be  described  (strictly 
speaking  only  in  the  limit  of  infinitely  long  chains)  using 
the  following  characteristic  ratio  Cn  [24] : 


Cn  = 


(8) 


where  /  is  the  bond  length  and  M—  1  the  number  of  bonds 
along  the  backbone  of  the  molecule. 

In  the  absence  of  short-range  interactions  between  differ¬ 
ent  segments  of  a  chain,  the  mean  square  end-to-end  length  of 
a  chain  is  defined,  in  accordance  to  the  random- wall  process, 
by  the  denominator  on  the  left  hand  side  in  Eq.  (8).  Hence, 
Cn  represents  a  ratio  of  the  actual  and  the  random- wall  based 
mean  square  end-to-end  chain  lengths. 


The  sum  on  the  right  side  of  Eq.  (9)  divided  by  Na  repre¬ 
sents  the  mean  square  displacement  (MSD)  of  species  a ,  Na 
is  the  number  of  diffusing  a  particles,  t  time  and  Ri(t)  is  the  ra¬ 
dial  displacement  of  particle  i  of  species  a  at  time  t.  It  must  be 
noted  that  Eq.  (9)  pertains  to  the  so-called  Einstein-diffusion 
and,  hence,  is  valid  only  when  diffusion  of  the  particles  is  a 
random  walk  process,  i.e.  the  displacement  of  a  particle  at 
some  time  is  uncorrelated  with  its  displacement  at  any  pre¬ 
vious  or  future  time.  In  this  regime  of  diffusion,  there  is  a 
linear  relation  between  the  mean  square  displacement  and 
the  time.  When  the  polymer-chain  surrounding  of  a  diffusing 
particle  inhibits  free  movement  of  such  particle,  i.e.  the  par¬ 
ticle  is  confined  (at  least  temporarily)  to  a  small  space,  the 
mass-transport  process  is  called  an  anomalous-diffusion.  In 
this  case  (\Ri(t)  —  7^-(0)|2)  a  tn,  where  n<  1,  and  Eq.  (9) 
is  not  applicable.  When  n  >  1  in  the  (|  Ri(t)  —  Ri( 0)|2)  a  tn 
relation,  on  the  other  hand,  the  mass-transport  is  dominated 
by  a  non-diffusion  process.  For  instance,  when  n-  2,  the  par¬ 
ticle  transport  takes  place  by  the  so-called  free-flight  process 
in  which  particles  move  in  a  low-density  medium  without 
any  interactions  with  the  medium.  To  identify  the  portion 
of  the  simulation  runs  within  which  the  Einstein-diffusion  is 
the  dominant  mass-transport  mechanism,  a  log(MSD)  versus 
log(0  plot  should  be  generated  and  the  slope  of  this  curve 
compared  with  unity. 

When  the  Einstein-diffusion  regime  is  reached,  the  ionic 
conductivity,  cr,  can  be  determined  using  the  following  Ein¬ 
stein  equation: 


6  t(V)kT 


Z?{[Ri(t)  -  Ri(0)f) 


+  2J2  ZiZj([Ri(t )  -  Rj(())][Rj(t)  -  Rj( 0)])  (10) 


]>l 


2. 7.  Transport  properties  analysis 

When  the  mass  or  the  ionic  transport  properties  of  the 
solid  polymer  electrolytes  are  considered  in  membrane  ap¬ 
plications,  three  fundamental  phenomena  must  be  analyzed: 

(1)  adsorption  of  the  transported  species  onto  the  membrane; 

(2)  their  transport  through  the  membrane;  (3)  desorption  of 


where  e  is  the  electronic  charge,  Z  the  valence,  k  the  Boltz¬ 
mann  constant  and  T  is  the  temperature.  The  first  term  on  the 
right-hand  side  of  Eq.  (10)  represents  the  sum  over  all  ionic 
mean  square  displacements  weighted  by  the  ionic  charges 
while  the  second  term  is  the  sum  of  correlations  of  the  ionic 
displacements  and  describes  the  contribution  of  the  electro¬ 
static  interactions  to  ionic  conductivity. 
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2.8.  Strength  analysis 


3.  Results  and  discussion 


When  the  PEO  sulfonic  acid-based  solid  electrolytes  are 
used  in  thin  membrane  applications,  they  must  generally 
satisfy  not  only  high  ionic-conductivity  requirements,  but 
also  meet  strength  requirements  in  order  to  ensure  a  long 
life  cycle.  To  assess  the  level  of  strength  in  the  PEO  sulfonic 
acid-based  materials  analyzed  in  the  present  work,  classi¬ 
cal  molecular  dynamics  calculations  are  employed  using  the 
Accelrys’  program  discover  [23].  The  molecular  dynamics 
procedure  utilized  involves  the  following  steps: 


(a)  The  atomic/molecular  structure  in  the  initial  (cubic,  un¬ 
deformed)  cell  is  first  equilibrated  by  treating  the  parti¬ 
cles  in  the  cell  and  the  cell  itself  as  an  NPT  ensemble, 
where  N  is  the  number  of  particles,  P  the  pressure  and 
T  the  absolute  temperature.  During  this  procedure,  the 
initial  particle  positions  and  the  (cubic)  cell  edge  length 
are  adjusted  to  keep  the  pressure  constant  and  equal  to  its 
target  value  (the  atmospheric  pressure  of  0.1  MPa  in  the 
present  case).  The  Anderson  method  for  pressure  control 
[23]  which  allows  a  change  in  the  cell  volume  but  not  in 
its  shape  is  used  in  this  step  of  the  molecular  dynamics 
simulations. 

(b)  NaT  molecular  dynamics  simulations  (where  a  is  the 
stress  tensor)  are  next  utilized  within  which  tensile  (nor¬ 
mal)  stress  is  applied  in  one  direction  at  a  constant  rate  of 
1  bar/ps  while  the  remaining  two  normal  stresses  are  kept 
at  their  0.1  MPa  level.  In  addition,  all  shear  stresses  are 
kept  at  a  zero  level.  The  magnitude  of  the  resulting  longi¬ 
tudinal  (normal)  strain  is  then  monitored  during  loading 
and  used  to  construct  the  tensile  stress- strain  curve. 

(c)  The  procedure  in  (b)  is  continued  until  “plastic”  yielding 
associated  with  an  abrupt  increase  in  the  longitudinal 
strain  rate  is  observed. 


An  example  of  the  deformed  cell  of  type  A  is  given  in 
Fig.  2(b). 


3.1.  Vibrational  spectral  analysis 

To  validate  the  NJPCFF  forcefield  used  in  the  present 
work,  the  local  mode  method  discussed  in  Section  2.5  is  ap¬ 
plied  to  an  isolated  PEO  sulfonic  acid  dianion  in  vacuum  and 
in  water.  The  results  obtained  are  compared  with  their  avail¬ 
able  corresponding  experimental  counterparts  for  the  PEO 
sulfonic  acid  (as  originally  compiled  by  Ennari  et  al.  [10-17]) 
and  a  summary  of  this  comparison  is  given  in  Table  1 .  Based 
on  the  results  displayed  in  Table  1,  the  following  main  find¬ 
ings  can  be  established: 

(a)  The  computed  vibrational  frequencies  for  the  sulfonic 
acid  group  in  PEO  sulfonic  acid  dianion  in  vacuum 
(1020  cm-1)  and  in  water  (1100  cm-1)  listed  in  Row 
10  are  quite  close  to  their  experimental  counterparts 
(1040-1090  cm-1),  considering  the  fact  that  the  accuracy 
of  the  local  mode  method  is  about  100  cm-1.  This  find¬ 
ing  suggests  that  the  additional  bonding  and  non-bonding 
forcefield  parameters  added  to  the  NJPCFF  to  account 
for  the  polymer/ion  interactions  correctly  account  for  the 
atomic-level  dynamics  of  the  sulfonic  end  groups  in  the 
PEO  sulfonic  acid  dianions. 

(b)  The  computed  and  measured  frequencies  cited  in  (a) 
are  clearly  different  than  the  measured  (850-910  cm-1) 
frequency  for  the  S— O  bond  stretch  (Row  8),  and  are 
also  different  than  the  measured  (1160  cm-1)  frequency 
for  the  S=0  bond  stretch  in  the  PEO  sulfonic  acid, 
Row  7. 

(c)  Since  the  difference  in  the  S— O-  bond  stretch  frequen¬ 
cies  in  the  PEO  sulfonic  acid  in  vacuum  and  in  water 
(Row  10)  are  within  the  range  of  accuracy  for  the  lo¬ 
cal  mode  method,  no  definite  conclusion  can  be  drawn 
regarding  the  effect  of  the  surrounding  medium  on  the 
atomic-level  dynamics  of  the  sulfonic  acid  end  group  in 
the  PEO  sulfonic  acid  dianions. 

(d)  A  comparison  of  the  remaining  computed  and  experi¬ 
mental  frequencies  shown  in  Table  1  indicates  that,  in 
general,  NJPCFF  predicts  quite  accurately  the  atomic- 


Table  1 

Computed  vibrational  absorption  frequencies  in  cm-1  for  the  PEO  sulfonic  acid  dianion  in  vacuum  and  water,  and  computed  and  measured  [3,16-20]  vibrational 
frequencies  for  PEO  sulfonic  acid 


Row 

PEO  sulfonic  acid  dianion 

In  vacuum  (computed) 

In  water  (computed) 

PEO  sulfonic  acid  (measured) 

Dominant  bond  assignment 

1 

3670 

3710 

O— H  in  water,  intra-molecular 

2 

3500 

3520 

3400-3600 

O— H  in  H30+,  hydrogen  bond 

3 

3550 

3400-3600 

O— H  in  PEO/water 

4 

2930 

2850 

2925 

C-H 

5 

1250 

1270 

800-1350 

c-c 

6 

1230 

1250 

1100 

c-o 

7 

1120-1170 

S—O 

8 

850-910 

s-o 

9 

1100 

600-700 

so3~  h3o+ 

10 

1020 

1100 

1070-1140 

s-cr 

11 

1090 

1090 

1070-1100 

c-s 
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Fig.  3 .  Radial  dependence  of  the  pair  correlation  function  between  the  proton 
(H+)  and  the  sulfur  atom  (S)  in  the  PEO  sulfonic  acid  dianion  and  of  the 
corresponding  coordination  number  for  the  same  type  of  particles. 

level  dynamics  of  the  PEO  sulfonic  acid  and  of  the  PEO 
sulfonic  acid  dianions. 


3.2.  Pair  correlation  and  coordination  number  analyses 

The  coordination  between  different  types  of  particles  is 
studied  by  analyzing  the  corresponding  pair  correlation  func¬ 
tions  g(r)  and  the  associated  coordination  numbers  n{r)  in  the 
first  and,  when  available,  in  the  second  solvation  shells.  An 
example  of  the  pair  correlation  function  versus  the  particles 
separation  distance  plot  is  given  in  Fig.  3.  A  summary  of  the 
results  of  the  coordination  analysis  is  given  in  Table  2.  These 
results  can  be  interpreted  as  follows: 

(a)  The  proton-  •  -ether-oxygen  (oxygen  in  the  backbone  of 
the  PEO  sulfonic  acid  and  the  PEO  polymer  chains)  cor¬ 
relation  function  and  the  single  solvation- shell  coordina¬ 
tion  number  decrease  slightly  while  the  corresponding 
solvation- shell  radius  increases  slightly,  from  cell  type 
A  to  cell  type  C,  Row  1.  This  finding  can  be  related  to 
the  number  of  water  molecules  in  the  cell  and  suggests 


that  the  H+  and  ether-0  particles  tend  to  become  more 
separated  and  less  coordinated  by  the  introduction  of  ad¬ 
ditional  water  molecules.  Due  to  the  large  dielectric  con¬ 
stant  of  water  and  the  observed  increased  H+-  •  -ether-0 
separation  distances  in  the  B  and  C  types  of  cells,  the 
non-bonding  (Coulomb  and  van  der  Waals)  interactions 
are  expected  to  decrease  when  the  material  at  hand  is 
hydrated. 

(b)  An  increased  number  of  water  molecules  in  the  B  and 
C  types  of  cells  also  appears  to  reduce  slightly  the  co¬ 
ordination  number  of  protons  within  the  first  solvation- 
shell  of  sulfur  atoms  in  the  sulfonic  acid  end  groups  of 
the  PEO  sulfonic  acid  dianions  and  to  increase  slightly 
the  solvation- shell  radius,  Row  2.  This  effect  is,  how¬ 
ever,  much  less  pronounced  in  the  case  of  the  second 
solvation- shell. 

(c)  Two  well-defined  solvation  shells  of  water  molecules 
(denoted  by  their  oxygen  atom,  02)  around  protons  with 
comparable  respective  radii  and  comparable  respective 
coordination  numbers  in  the  three  types  of  cells  are  ob¬ 
served,  Row  3.  This  finding  suggests  that  protons  are 
mainly  found  in  water-rich  regions  located  at  relatively 
large  distances  from  the  PEO  sulfonic  acid  dianions  and 
the  PEO  polymer  molecular  chains.  Hence,  protons  are 
expected  to  be  quite  mobile  in  the  presence  of  water  in 
the  materials  at  hand,  and  consequently,  the  contribution 
of  the  hopping  mechanism  to  proton  transport  is  expected 
to  be  significant. 

(d)  Coordination  between  ether-oxygen  and  oxygen  (named 
03)  in  the  hydronium  ions,  Row  4,  is  seen  to  result  in 
the  formation  of  a  strong  second  solvation-shell  whose 
strength  and  the  radius  seem  to  increase  slightly  as  the 
number  of  water  molecules  increases  from  cell  type  A  to¬ 
ward  cell  type  C.  This  finding  suggests  that  ether-oxygen 
atoms  and  hydronium  ions  are  separated  by  water 
molecules,  and  that  the  separation  distance  increases  with 
the  number  of  water  molecules  present  in  the  system. 

(e)  An  increase  in  the  number  of  water  molecules  in  cells  B 
and  C  gives  rise  only  to  a  slight  increase  in  the  radius  of 
the  first  solvation- shell  for  hydronium  ions  surrounding 
the  sulfur  atom  in  the  sulfonic  acid  end  group,  Row 
5.  The  corresponding  coordination  number,  however, 


Table  2 


The  results  of  the  pair  correlation  and  the  coordination  number  analyses  carried  out  in  the  present  work 


Row  number 

Pair  correlation  function 

Cell  A 

Cell  B 

Cell  C 

Radius  (A) 

Coord,  number 

Radius  (A) 

Coord,  number 

Radius  (A) 

Coord,  number 

1 

H+-  • 

•ether-O 

2.44 

2.05 

2.51 

1.91 

2.55 

1.90 

2 

H+-  • 

•S 

3.54  (6.10) 

2.0  (2.32) 

3.59  (6.07) 

1.79  (2.40) 

3.6(6.11) 

1.75  (2.09) 

3 

H+-  • 

•02 

2.41  (4.91) 

2.95  (8.25) 

2.45  (5.25) 

2.66  (9.11) 

2.56  (4.99) 

2.67  (8.72) 

4 

03-  • 

•ether-O 

3.42  (4.52) 

0.41  (2.95) 

3.51  (4.69) 

0.40  (3.17) 

3.54  (4.81) 

0.35  (3.19) 

5 

03-  • 

•S 

4.13(6.31) 

1.21  (2.09) 

4.18(6.49) 

1.25  (2.00) 

4.19(6.58) 

1.12(1.95) 

6 

02-  • 

•03 

3.15 

7.95 

3.29 

8.51 

3.49 

8.49 

7 

02- • 

•ether-O 

3.42  (4.12) 

1.09  (2.12) 

3.49  (4.19) 

1.02(1.80) 

3.59  (4.30) 

1.01  (1.59) 

8 

02- • 

•S 

4.12(5.96) 

10.62  (16.32) 

4.22  (6.22) 

10.52  (15.99) 

4.63  (6.12) 

10.54  (16.92) 

The  number  within  parenthesis  pertain  to  the  second  solvation-shell.  02:  oxygen  atom  in  H2O;  03:  oxygen  atom  in  H30+. 
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while  being  affected  by  the  number  of  water  molecules 
in  the  cell,  varies  with  this  number  in  a  more  complex 
manner.  A  slightly  larger  effect  of  hydration  is  seen  in 
the  case  of  the  second  solvation- shell. 

(f)  Water  molecules  form  a  strong  single  solvation- shell 
around  the  hydronium  ions  and,  while  the  strength  of 
this  solvation-shell  clearly  varies  with  the  number  of 
water  molecules  in  the  cell,  Row  6,  the  nature  of  this 
interaction  appears  to  be  quite  complex. 

(g)  In  the  case  of  all  three  types  of  cells,  water  molecules 
are  strongly  correlated  with  the  ether-oxygen  atoms 
which  results  in  the  formation  of  two  well-defined  H2O 
solvation  shells,  Row  7.  The  radii  and  the  coordination 
numbers  of  these  shells,  however,  are  not  significantly 
affected  by  differences  in  the  number  of  water  molecules 
between  the  three  types  of  cells. 

(h)  Water  molecules  are  strongly  correlated  with  sulfur 
atoms  in  the  sulfonic  acid  end  groups  resulting  in  the 
formation  of  two  well-defined  solvation  shells  of  H2O 
molecules,  Row  8.  It  is  generally  believed  (e.g.  [27]) 
that  at  least  10  water  molecules  should  be  correlated 
with  each  sulfonic  acid  end  group  in  order  to  obtain 
good  ionic  conductivity  in  the  PEO  sulfonic  acid-based 
solid  polymer  electrolytes.  Based  on  the  water-sulfur 
correlation  data  listed  in  Row  8,  Table  2,  one  can  expect 
that  the  materials  associated  with  all  three  types  of  cells 
are  ionically  conductive. 

3.3.  Conformation  analysis 

The  structure  of  the  PEO  sulfonic  acid  dianion  and  PEO 
polymer  chains  is  studied  by  analyzing  the  distributions  of  the 
OSCC,  SCCC,  CCCO,  CCOC  and  OCCO  dihedral  angles. 
The  dihedral-angle  for  a  chain  of  three  linearly-connected 
bonds  where  the  bonds  are  represented  with  vectors  r;_i,  r; 
and  r/+i,  respectively,  is  the  angle  between  the  plane  of  the 
vectors  r;_  1  and  q,  and  the  plane  of  the  vectors  77  and  77+1 .  A 
typical  (SCCC)  dihedral- angle  distribution  plot  is  shown  in 
Fig.  4.  The  chain  conformations  associated  with  a  dihedral- 
angle  of  0°  are  referred  as  the  trans  (T)  conformation  states 
while  those  associated  with  dihedral  angles  of  —120°  and 
120°  are  denoted  as  the  gauche  minus  (G“)  and  the  gauche 
plus  (G+)  conformation  states,  respectively,  is  shown  in  Fig.  4. 

A  summary  of  the  dihedral- angle  results  obtained  in  the 
present  work  is  given  in  Table  3.  To  enable  a  simple  quan¬ 
titative  comparison  of  the  populations  of  the  three  dihedral- 
angle  conformation  states  (T,  G+  and  G-)  in  the  three  types 
of  cells,  the  T  conformation  state  is  taken  to  correspond  to 
a  —60°  to  60°  dihedral-angle  range,  the  G-  conformation 
state  to  a  —180°  to  60°  dihedral-angle  range  while  the  G+ 
conformation  state  to  a  60°-180°  dihedral-angle  range.  Due 
to  torsional  symmetry  of  the  PEO  sulfonic  acid  dianion  and 
the  PEO  molecules,  the  G_  and  G+  conformation  states  are 
generally  found  to  be  present  with  approximately  equal  prob¬ 
abilities  for  each  of  the  five  dihedral  angles  analyzed.  There¬ 
fore,  in  Table  3,  only  the  percentages  of  the  T  conformation 


Fig.  4.  Probability  density  plot  for  the  average  SCCC  dihedral-angle  for  the 
PEO  sulfonic  acid  dianion. 

state  for  the  five  dihedral  angles  are  given.  The  corresponding 
percentages  of  the  G_  and  the  G+  conformation  states  can  be 
readily  computed. 

The  results  displayed  in  Table  3  can  be  summarized  as 
follows: 

(a)  The  mean  values  for  the  corresponding  five  dihedral  an¬ 
gles  are  very  similar  in  the  three  types  of  cells  and,  thus, 
not  significantly  affected  by  the  existing  differences  be¬ 
tween  the  three  types  of  cells  in  the  amount  of  PEO  poly¬ 
mer  and  the  number  of  water  molecules. 

(b)  The  T,  G_  and  G+  conformation  states  of  the  OSCC 
dihedral-angle  are  almost  equally  populated  which  is 
most  likely  the  result  of  symmetry  of  the  PEO  sulfonic 
acid  dianion  end  groups. 

(c)  The  SCCC  dihedral-angle  is  almost  completely  populated 
with  the  trans  conformation  state. 

(d)  The  T,  G-  and  G+  conformation  states  for  the  CCCO 
dihedral-angle  are  nearly  equally  populated. 

(e)  The  conformation  of  the  CCOC  dihedral-angle  is  domi¬ 
nated  by  its  T  state  although  the  G-  and  the  G+  confor¬ 
mation  states  are  present  in  measurable  quantities. 

(f)  The  trans  state  of  the  OCCO  dihedral-angle  is  deficient 
(i.e.  its  fraction  is  smaller  than  33%)  indicating  the  occur¬ 
rence  of  the  so-called  gauche  effect.  This  effect  is  gener- 

Table  3 

Mean  values  of  the  percent  of  the  trans  conformation  state  in  various  dihedral 
angles  in  the  PEO  sulfonic  acid-based  polymeric  electrolyte  studied  in  the 
present  work 


Dihedral-angle 

Cell  A 

Cell  B 

Cell  C 

OSCC 

32.8 

34.1 

33.8 

SCCC 

94.1 

95.0 

93.3 

CCCO 

34.0 

32.1 

33.7 

CCOC 

82.2 

85.8 

84.7 

OCCO 

30.1 

27.4 

28.3 
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ally  promoted  by  high  polarity,  high  dielectric  constant 
and  the  presence  of  hydrogen  bonding  in  the  surrounding 
medium  and  all  these  conditions  are  present  in  the  mate¬ 
rial  at  hand  [14] .  In  polymers  with  a  pronounce  gauche  ef¬ 
fect,  a  new  cation  transport  mechanism  is  often  observed 
in  which  the  cations  hop  from  one  ether-oxygen  to  the 
adjacent  one.  However,  such  a  mechanism  is  not  expected 
to  play  a  major  role  in  material  analyzed  in  the  present 
work  considering  the  fact  that  no  significant  correlation  is 
found  between  the  cations  (protons  and  hydronium  ions) 
and  the  ether-oxygen  atoms  residing  in  the  backbone  of 
PEO  sulfonic  acid  and  the  backbone  of  PEO  molecules. 

It  should  be  noted  that  the  results  displayed  in  Table  3  re¬ 
veal  that  the  number  of  water  molecules  in  the  computational 
cell  does  not  have  a  major  effect  on  the  population  distribution 
of  various  dihedral  angles,  and  that  the  variations  observed 
are  a  reflection  of  the  statistical  nature  of  the  problem  at  hand. 

3.4.  Dimensional  analysis 

To  get  a  further  insight  into  the  conformation  of  the  PEO 
sulfonic  acid  dianions  in  the  three  cells,  Eqs.  (6)— (8)  are  used 
to  compute  the  radius  of  gyration,  the  end-to-end  distance  and 
the  characteristic  ratio  for  the  PEO  sulfonic  acid  dianions  in 
the  three  types  of  cells.  The  results  of  this  calculation  are 
given  in  Table  4  and  can  be  summarized  as  follows: 

(a)  The  mean  values  of  the  corresponding  dimensional  pa¬ 
rameters  do  not  differ  significantly  in  the  three  types  of 
cells. 

(b)  The  mean  values  for  the  radius  of  gyration  and  for  the 
end-to-end  distance  in  all  three  types  of  cells  are  smaller 
by  a  factor  of  approximately  60%  than  their  counterparts 
for  PEO  sulfonic  acid  dianions  in  vacuum  [14] .  This  find¬ 
ing  suggests  that  these  dianions  are  more  curved  or  coiled 
in  the  three  cells  studied  in  the  present  work  and  more 
extended  when  present  as  isolated  molecules  in  vacuum. 
At  least  two  reasons  may  be  cited  for  the  observed  change 
in  the  PEO  sulfonic  acid  dianion  conformations:  (i)  due 
to  a  high  value  of  the  dielectric  constant  of  water,  the  re¬ 
pulsion  between  the  two  negatively  charged  ends  of  the 
PEO  sulfonic  acid  dianions  is  reduced  in  the  case  of  the 
three  types  of  cells  analyzed  in  the  present  work.  This 
enables  the  dianions  ends  to  get  closer  to  each  other;  and 
(ii)  the  observed  gauche  effect  results  naturally  in  the 
chains  with  lower  values  of  their  radius  of  gyration  and 
of  their  end-to-end  distance. 

Table  4 

Average  values  and  standard  deviations  for  the  radius  of  gyration,  the  end-to- 
end  distance  and  the  characteristic  ratio  for  the  PEO  sulfonic  acid  dianions 
in  the  three  cells  studied  in  the  present  work 


Dimensional  parameter 

Cell  A 

Cell  B 

Cell  C 

Radius  of  gyration  (A) 

6.22  ±  0.14 

5.98  ±  0.22 

6.11  ±  0.25 

End-to-end  distance  (A) 

12.3  ±  3.4 

11.9  ±  4.4 

12.5  ±  3.3 

Characteristic  ratio  (N/A) 

3.2  ±  0.1 

3.3  ±  0.1 

3.0  ±  0.1 

Table  5 

The  average  values  (in  1 0-9  m2/s)  and  the  standard  deviations  for  the  diffu¬ 
sion  coefficient  of  various  species  in  the  PEO  based  hydrated  solid  electrolyte 
studied  in  the  present  work 


Diffusing  species  Cell  A  Cell  B  Cell  C 


H+ 

0.36 

± 

0.04 

0.46 

± 

0.05 

0.53 

± 

0.04 

h3o+ 

0.39 

± 

0.04 

0.45 

± 

0.06 

0.50 

± 

0.07 

h2o 

0.82 

± 

0.06 

0.86 

± 

0.07 

0.92 

± 

0.10 

CT 

0.03 

± 

0.02 

0.04 

± 

0.02 

0.04 

± 

0.01 

Ether-O 

0.02 

± 

0.01 

0.02 

± 

0.01 

0.02 

± 

0.01 

(c)  The  values  for  the  characteristic  ratio  listed  in  Table  4 
are  also  lower  than  their  counterparts  for  the  isolated 
PEO  sulfonic  acid  dianions  in  vacuum  ( Cn  =  9A ,  [14]) 
confirming  the  observation  made  above  that  hydration 
makes  these  ions  less  extended. 

It  should  be  noted  that  the  results  presented  in  Table  4 
suggest  that  the  radius  of  gyration,  the  end-to-end  distance 
and  the  characteristic  ratio  are  not  significantly  affected  by 
the  number  of  water  molecules  in  the  cell,  as  evidenced  by 
relatively  small  changes  in  the  magnitudes  of  these  quantities 
due  to  an  increase  in  the  number  of  water  molecules  compared 
with  the  standard  deviation  values  at  a  given  number  of  water 
molecules. 

3.5.  Mass  and  charge  transport  properties 

Diffusion  coefficients  for  the  proton,  the  hydronium  ion, 
the  water,  the  oxygen  in  the  sulfonic  acid  end  group  and  the 
ether-oxygen  of  the  PEO  sulfonic  acid  and  PEO  polymer 
chains  for  the  three  types  of  cells  computed  using  Eq.  (9) 
are  summarized  in  Table  5.  The  condition  for  the  Einstein- 
diffusion  discussed  in  Section  2.7  was  met  in  the  case  of  all 
particles  in  all  three  types  of  cells.  The  results  displayed  in 
Table  5  can  be  briefly  summarized  as  follows: 

(a)  The  diffusion  coefficients  of  the  proton,  the  hydronium 
ion  and  the  water  are  mutually  comparable  in  the  case 
of  each  of  the  three  types  of  cells.  The  diffusion  coeffi¬ 
cients  of  the  oxygen  ion  in  the  sulfonic  end  group  and 
of  the  ether-oxygen  in  the  backbone  of  the  PEO  sul¬ 
fonic  acid  and  the  PEO  polymer  chains  are  smaller  by 
at  least  an  order  of  magnitude.  This  finding  can  be  read¬ 
ily  understood  considering  the  fact  that  H+,  H30+  and 
H2O  are  small  free  particles  while  O-  and  ether-0  are 
covalently-bonded  constituents  of  the  PEO  sulfonic  acid 
dianion  chains  while  ether-0  is  also  a  constituent  of  the 
PEO  molecular  chains. 

(b)  In  all  three  types  of  cells,  the  diffusion  coefficient  of  wa¬ 
ter  is  the  highest  and  increases  slightly  with  the  number 
of  water  molecules  in  the  cell.  While  this  finding  can 
be  partly  attributed  to  a  relatively  large  (40  wt.%)  weight 
percent  of  H2O  in  the  system,  it  may  also  suggest  that  the 
other  two  small  particles,  i.e.  H+  and  H30+  cations,  are 
involved  in  non-bonding  interactions  with  the  PEO  sul- 
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fonic  acid  dianions  and  the  PEO  molecules  which  lowers 
their  diffusivities. 

(c)  As  the  number  of  water  molecules  increases  from  the  cell 
type  A  to  the  cell  type  C,  the  diffusion  coefficients  of  H+, 
H30+  and  H2O  all  increase.  Between  the  two  cations,  the 
diffusion  coefficient  of  H+  increases  at  a  somewhat  faster 
rate  than  that  of  H30+  as  the  number  of  water  molecules 
is  increased.  This  finding  suggests  that  as  the  number  of 
water  molecules  increases,  the  role  of  the  proton  hopping 
mechanism  in  the  H+  transport  and,  hence,  the  magnitude 
of  the  ionic  conductivity  increase. 

(d)  Since  the  diffusion  coefficients  of  O-  and  of  ether-0 
are  substantially  lower  than  those  of  H+  and  H30+,  the 
latter  particles  appear,  for  the  most  part  to  be  free  to 
move  through  the  material  as  opposed  to  being  confined 
by  or  bonded  to  and  moving  with  the  PEO  sulfonic  acid 
dianion  and  the  PEO  molecular  chains. 

The  total  ionic  conductivity  computed  using  Eq.  (10)  are: 
6.4,  6.9  and  7.2mS/cm  for  the  cell  types  A,  B  and  C,  re¬ 
spectively.  The  standard  deviation  for  these  results  is  around 
0.8mS/cm.  For  comparison,  the  corresponding  experimen¬ 
tal  values  obtained  using  impedance  measurements  are  1.5 
and  0.015  mS/cm  at  the  relative  humidity  levels  of  75.3  and 
38.2%,  respectively  [8].  While  the  information  provided  in 
Ref.  [8]  does  not  enable  determination  of  the  correspond¬ 
ing  water  weight  percents  in  the  solid  polymer  electrolyte, 
the  computed  values  appear  to  be  somewhat  higher  than  ex¬ 
pected.  One  possible  explanation  for  this  shortcoming  of  the 
model  is  the  way  in  which  the  proton  hopping  mechanism 
was  modeled.  That  is,  proton  hopping  between  adjacent  wa¬ 
ter  molecules  is  assumed  to  take  place  without  breaking  of 
the  chemical  bonds,  the  process  which  is  generally  associ¬ 
ated  with  a  large  activation  energy  and,  hence,  a  lower  rate. 
Consequently,  the  computed  contribution  of  proton  hopping 
to  the  ionic  conductivity  in  the  material  at  hand  is  expected 
to  be  somewhat  overestimated. 

By  analyzing  separately  the  contributions  of  the  protons 
and  the  hydronium  ions  to  the  ionic  conductivity  of  the  mate¬ 
rial  at  hand,  it  is  established  that  in  all  three  types  of  cells  these 
two  contributions  are  comparable.  However,  as  the  number 
of  water  molecules  is  increased  from  Cell  A  toward  Cell  C, 
the  relative  contribution  of  the  proton  hopping  mechanism  to 
the  ionic  conductivity  increases. 

3.6.  Mechanical  strength  analysis 

The  average  tensile  stress-strain  curves  along  with  the  cor¬ 
responding  one-standard  deviation  error  bars  generated  using 
the  procedure  described  in  Section  2.8  for  the  three  types  of 
cells  are  shown  in  Fig.  5.  The  corresponding  experimental 
stress-strain  curve  for  a  40  wt.%  hydrated  PEO  polymer  [28] 
is  also  shown  in  Fig.  5  for  comparison.  The  results  displayed 
in  Fig.  5  can  be  summarized  as  follows: 

(a)  Despite  a  considerable  scatter  in  the  simulation  results, 
the  computed  tensile  stress-strain  curves  for  the  materials 


Fig.  5 .  Experimental  stress-strain  curve  for  a  40  wt.  %  hydrated  PEO  polymer 
[28]  and  computed  stress-strain  curve  for  the  three  types  of  cells  containing 
PEO  sulfonic  acid  solid  polymer  electrolyte. 


associated  with  the  three  types  of  cells  are  surprisingly 
quite  close  to  the  experimental  curve. 

(b)  In  spite  of  a  significant  value  of  the  standard  deviation  for 
the  longitudinal  strain  in  the  three  types  of  cells,  material 
associated  with  the  cell  type  B  is  found  to  possess  statis¬ 
tically  the  largest  value  of  its  yield  strength.  This  finding 
suggests  that  at  a  same  level  of  hydration,  there  may  be 
an  optimum  fraction  of  the  PEO  polymer  which  gives 
the  best  combination  of  ionic  conductivity  and  strength 
in  the  PEO  sulfonic  acid-based  solid  polymer  electrolyte. 
This  conjunction  will  be  explored  in  more  details  in  the 
future  work. 

(c)  A  detailed  examination  of  the  conformations  of  the  PEO 
sulfonic  acid  dianions  and  of  the  PEO  polymer  molecules 
reveals  that  numerous  gauche- to- trans  conformation- 
state  transitions  in  the  OCCO  dihedral- angle  take  place 
during  yielding.  Consequently,  the  mean  values  of  the 
percent  of  the  trans  conformation  state  in  the  three  types 
of  cells  at  a  strain  level  of  10%  are  found  to  be  in  a 
39.6-43.2%  range.  This  values  are  substantially  higher 
than  their  counterpart  (27.4-30.1%,  Table  3)  in  the  ini¬ 
tial  state.  The  gauche-to-tra^s  conformation-state  transi¬ 
tions,  thus,  appear  to  be  the  primary  yielding  mechanism 
in  the  material  under  investigation. 

4.  Conclusions 

Based  on  the  results  obtained  in  the  present  work,  the 

following  main  conclusions  can  be  drawn: 

1 .  Molecular  modeling,  in  conjunction  with  experimental  in¬ 
vestigations,  can  be  a  very  instrumental  tool  in  identifying 
the  microstructural  features  and  the  atomic-level  mech- 
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anisms  which  control  various  physical  and  mechanical 
properties  of  solid  polymer  based  electrolytes. 

2.  In  PEO  sulfonic  acid  dianion  based  solid  polymer  elec¬ 
trolytes,  both  the  classical  cation  diffusion  and  a  proton 
hopping  mechanism  play  important  roles  in  the  ionic  con¬ 
duction. 

3.  As  expected,  hydration  promotes  ionic  conductivity  but 
degrades  the  strength  of  the  PEO  sulfonic  acid  solid  elec¬ 
trolytes. 

4.  At  a  fixed  level  of  hydration,  it  appears  that  there  is  an 
optimum  level  of  the  PEO  polymer  content  which  gives 
the  best  combination  of  the  material  ionic  conductivity 
and  strength. 
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